www.gusucode.com > matlab编程NSCT分解 图像融合 各个融合指标评价体系 分解源码程序 > matlab编程NSCT分解 图像融合 各个融合指标评价体系 分解源码程序/NSCT/nsct1.m

    function  coeffs=nsct(im,n)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Image decomposition by nonsubsampled contourlet transform (NSSC).
% This is the iterated filter bank that computes the nonsubsampled
% contourlet transform.  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% Parameteters:
nlevels = n ;        % Decomposition level
%pfilter = 'pkva' ;              % Pyramidal filter
dfilter = 'dmaxflat7'; %'cd' ;              % Directional filter
pfiltr= '9-7';
% Nonsubsampled Contourlet decomposition
%im=imread('zoneplate.png');
coeffs = nsctdec( double(im), nlevels ,dfilter,pfiltr);


% disp( nlevels); disp(dfilter);disp(pfiltr);
% 
% % Display the coefficients
% %disp('Displaying the contourlet coefficients...') ;
% %shownsct( coeffs ) ;
% 
% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% % Nonsubsampled Contourlet transform (NSSC) reconstruction.
% % This is the inverse of nsscdec, i.e.
% % imrec = nsscrec(coeffs, dfilter, pfilter);
% % would reconstruct imrec = im
% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% 
% % Reconstruct image
% imrec = nsctrec( coeffs, dfilter, pfiltr ) ;
% 
% 
% % Show the reconstruction image and the original image
% %if 0
% figure;
% subplot(1,2,1), imagesc( im, [0, 255] ); 
% title('Original image' ) ;
% colormap(gray);
% axis image off;
% subplot(1,2,2), imagesc( imrec, [0, 255] );
% title('Reconstructed image' ) ;
% colormap(gray);
% axis image off;
% %end
% 
% mse = sum( sum( (imrec - double(im)).^2 ) );
% mse = mse / prod(size(im));
% 
% disp( sprintf('The mean square error is: %f', mse ) );
% disp(' ');